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Q ABSTRACT 

We present a rest-frame spectral stacking analysis of ~ 1000 X-ray sources detected in the XMM-COSMOS field in 
^ order to investigate the iron K line properties of active galaxies beyond redshift 2 ~ 1. In Type 1 AGN that have a 

typical X-ray luminosity of Lx ~ 1.5 x lO"*" (crgs^^) and z ~ 1.6, the cold Pe K at 6.4 keV is weak {EW ~ 0.05 keV), 
I in agreement with the known trend. In contrast, high-ionization lines of Fe xxv and Pc xxvi arc pronounced. These 

O high-ionization Fe K lines appear to have a connection with high accretion rates. While no broad Fc omission is detected 

in the total spectrum, it might be present, albeit at low significance (~ 2a), when the X-ray luminosity is restricted 
to the range below 3 x 10**^ ergs ^, or when an intermediate range of Eddington ratio around A ~ 0.1 is selected. 
I ^ , In Type II AGN, both cold and high-ionzation lines become weak with increasing X-ray luminosity. However, strong 

high-ionization Pe K {EW ~ 0.3 keV) is detected in the spectrum of objects at z > 2, while no 6.4 keV line is found. 
I It is then found that the primary source of the high-ionization Pe K emission is those objects detected with Spitzer- 

^ MIPS at 24fj,m. Given their median redshift of 2 ~ 2.5, their bolometric luminosity is likely to reach lO^^L© and the 

— . MlPS-detected emission most likely originates from hot dust heated by embedded AGN, probably accreting at high 

Eddington ratio. These properties match those of rapidly growing black holes in ultra-luminous infrared galaxies at the 
interesting epoch (2 ~ 2-3) of galaxy formation. 

(N Key words. X-rays: galaxies - Galaxies: active - Surveys 

1. Introduction exposure observations are available, e.g., Comastri et al 

(2011), Feruglio et al (2011), Norman et al (2002). 

>The iron K-shell (Fe K) emission at 6-7 keV is the most ttt-.i .i 1 i- 1 1 rj.i 1 / r, j 2\ 
, . , . 1 1- • xi V i. r 1 With the relatively large area or the sky (~2-aeg ) cov- 

prominent spectral line m the X-ray spectra or active galax- j 1 ,1 r^r^mirAa /o -n ^ 1 r,r^r^%\ -^1 

K> . J ^ ir - u-i ered by the COSMOS survey (Scoville et al 2007) with 

\/\ les and serves as a useful probe of various physical proper- -^^ , , , , ^ s ■ ^ a i i 1 viv m r at 

,. r 1 i- 1 • TT • J.T J. -1 ui its depth down to lamt nuxcs probed by XMM-Ncwton 

53 ties ol galactic nuclei. However, given the current available „ ,. , 1 r,r,nn\ • ^ i n • v 

C^.^°,^. V , '? r ^- 1,. (Cappelluti et al 2009), a sigmfacant number 01 luminous X- 

mstrumcntation, X-ray observations 01 active galactic mi- ^ , ■, 1 r ^ ■ ^ a ^-nt 1 ; ; 1 

1 • r\r^-\j\ ■ ii 1^1-^1 • • ^ ray sources at > 1, most 01 which are AGN, are detected, 

clci AGIN are still m the photon-limited regime m most ~r ^ i- • • i- • 1 i • t 

^ , , , J r n 1- - I- J i 1 While detection 01 a le K line m individual sources is lim- 

cases, and the study ot le K lines is limited to nearby ., , , n r rj.iv i • ^ /at • • • 

1 or fw r .1 1 1 ited to a very small traction ot the X-ray obiects IVlaimeri 

bright objects. For example, - 85 % of the unobscured 2007), we have made use of the large XMM-COSMOS 

AGN sample ol Bianchi et al (GAIXA, 2009), which is , , , .-.i . , , onn-rN j f i j. i • i. u 

, , ,. r 1 J Ar^AT • 4-1 vnTn/r dataset (Hasmger et al 2007) and spectral stacking to ob- 

the largest compilation ol unobscured AGN m the XMM- . ■ ■ . - i c c ■ ■ ■ c j 

. , . ^ , ji-r, ^ r^ A r> J i taiu integrated signals ol the taint Fe emission tor under- 

Newton archive, he at redshift z < 0.4. Beyond 2 = 1, . j- .i i • i r .• i • . i • i 

,. , . r J.- T • • i 11 standing the physical properties oi active galaxies at high 

the le K line properties ot active galaxies are virtually un- j.g(jgiiif^ r- j t- o o 

explored apart from a few exceptions for which very long 

Here we first summarize a few relevant facts concerning 

the Fc K line, obtained from observations of nearby AGN: 

* Email: kazushi.iwasawa@icc.ub.edu 1) The emission feature is dominated by the 6.4 keV line (up 
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to Fe xvii), indicating that a cold medium illuminated by 
an AGN is the primary source of the line emission. High- 
ionization lines from Fe XXV at 6.7 keV and Fe xxvi at 
6.97 kcV arc also seen in good quality data (e.g., Bianchi 
et al 2007) but arc much weaker than the cold line; 2) In 
unobscurcd AGN, a line profile with a redward extension is 
observed in a significant fraction of nearby Seyfert galaxies 
(de la Calle Perez et al 2010; Pukazawa et al 2011; Bianchi 
et al 2009; Nandra ct al 2007). Its origin has been under 
debate, but it may be a result of spectral distortion due 
to strong gravity, and therefore could serve as a powerful 
probe of the accretion flow in the relativistic region in the 
vicinity of a black hole (e.g., Fabian et al 2000); and 3) 
In heavily obscured AGN, observed X-rays arc suppressed, 
and thus reprocessed light is relatively enhanced. This leads 
to a strong appearance of Fe K in an observed spectrum 
(Awaki ct al 1991; Ghisellini, Haardt & Matt 1994; KroHk, 
Madau & Zycki 1994). A large equivalent width {EW ^ 1 
keV) of Fc K is then considered to be a characteristic feature 
of heavily obscured AGN with an absorbing Thomson depth 
of unity or larger, i.e., Compton-thick AGN. 

The goal of this work is to investigate more power- 
ful AGN at higher redshifts, as detected in the XMM- 
COSMOS field, to see whether the above properties hold or 
evolve with increasing luminosity /redshift. We show below 
that the Fe K properties of AGN in the XMM-COSMOS 
field appear different from those in nearby objects, and dis- 
cuss possible drivers of these changes. 

After presenting an overview of the XMM-COSMOS 
dataset and the spectral stacking method employed in this 
paper, results are presented for the whole dataset as well 
as its several subsets in order to investigate the variations 
of the Fe K line properties. 

The cosmology adopted here is Hq = 70 km s~^ Mpc~^, 

= 0.72, Om = 0.28. 

2. Observation and data reduction 

About 1800 point-like X-ray sources were detected in the 
XMM-COSMOS survey (Cappelluti et al 2009). The mean 
exposure time over the survey area is about 40 ks. Their op- 
tical identifications and redshifts compiled from literature, 
various spectroscopic campaigns (Lilly et al 2009; Trump 
et al 2009) and photometric redshifts (Salvato et al 2009, 
2011) are presented in Brusa et al (2010). The spectrum 
extraction procedure for each XMM-COSMOS source is de- 
scribed in Mainieri et al (2007), and we use those spectra 
obtained from the EPIC pn camera. 

We use the X-ray data analysis software packages 
HEASARC's FTOOLS 6.1 and XSPEC version 12 for the 
analysis presented below. 

3. The Sample 

The X-ray sources used in our spectral stacking analysis 
are selected from the XMM-COSMOS catalogue of Brusa 
ct al (2010). 1449 cxtragalactic sources with a secure optical 
counterpart ("flagID = 1" in Brusa et al 2010), and either a 
spectroscopic or photometric redshift, were selected for fur- 
ther filtering based on the quality of their X-ray data. Given 
our interest in the rest-frame 2-10 keV band, background 
corrected counts of each source in the rest frame 2-10 kcV 
band are used to discard poor quality spectra, which would 



only increase the noise in a final stacked spectrum. The dis- 
tribution of the counts is shown in Fig. 1. A clear break is 
seen at 20 counts, which would not appear if observed with 
infinite sensitivity and thus this break can be identified as 
the sensitivity limit (more correctly, below this limit, the 
uncertainty in background subtraction becomes too large, 
and returns in data stacking diminish). We adopted this 
break as a threshold for discarding poor quality spectra of 
406 sources. On the other end, we identified 23 sources with 
more than 400 counts as bright outliers (see Fig. 1), and 
we have separated them from the general stacking analysis. 
They are hereafter referred to as the "Bright 23" sample 
(see § 3.3). The exclusion of these brightest objects leaves 
1020 som-ces with 20-400 source counts. We refer to these 
1020 sources as the main sample, and the results shown be- 
low are derived from this sample unless stated otherwise. 
The total counts in the rest-frame 2-10 keV for all sources 
in the main sample is 77039 cts, half of which come from 
sources with less than 100 counts individually. For sources 
in the main sample, the typical source fraction of the accu- 
mulated rest-frame 2-10 keV counts is 50%, i.e., source 
and background counts are comparable. 

3.1. Main sample 

In the main sample, 643 sources have spectroscopic reel- 
shifts, and photometric redshifts arc available for the other 
377 sources. The median values of redshift are itot = 1-30 
for all sources, Zgp = 1.17 for the sources with spectroscopic 
redshift and Zph = 1.59 for the sources with photometric 
redshift. 

For sources with spectroscopic redshifts, those that are 
broad-line AGN based on their optical spectra, are assigned 
to "Class I" , while the remainder are assigned to "Class 
11". Salvato et al (2009, 2011) give the best-fit SED type for 
individual sources with photometric redshifts. The photo- 
metric sample is divided into "SED I" when the SED class 
is 19 or larger (those of broad-linc AGN, see Salvato et al 
2009 for details), and "SED 11" otherwise. 

For the spectral stacking discussion below, we divide 
the sample into two "Types" that are defined as follows: 
Type I = Class I + SED I, and Type II = Class II + 
SED II. In other words, Type I objects are broad-line AGN 
or objects with a SED similar to that of broad-line AGN, 
while Type II objects are all of the rest, including AGN with 
Seyfert 2 excitation, emission-line galaxies, and absorption- 
line galaxies. The Type I sample contains 508 sources with 
a median redshift of z = 1.60, while the Type II sample 
contains 512 sources with median redshift of 2; = 0.97. Their 
redshift distributions are shown in Fig. 2. Plots of the 2-10 
keV luminosity (with no correction for absorption) versus 
redshift for Type I and Type II objects in the main sample 
as well as the Bright 23 sample are shown in Fig. 3. 

Fig. 4 shows the 2-10 keV luminosity distributions for 
the Type I and Type II objects. The median values of log- 
arithmic luminosities are 44.16 and 43.65 (erg s~^), respec- 
tively. 

3.2. Broad-line AGN: Mbh and A 

For spectroscopically confirmed broad-line AGN (Class I), 

their black hole masses estimated from the optical broad 
lines are available from the literature (Trump et al 2009; 
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Table 1. Properties of the main and bright samples. 
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Spectro-z 


643 


1.17 


43.87 


0.71 ± 0.03 


57538 




O t 1 


1.59 


43.53 


0.34 ±0.05 




Type I 


508 


1.60 


44.16 


0.80 ±0.05 


48218 


Type II 


512 


0.97 


43.65 


0.21 ±0.05 


28821 






Type I objects 






Class I 


366 


1.57 


44.21 


0.80 ±0.03 


40022 


SED I 


142 


1.78 


44.03 


0.59 ±0.07 


8196 






Type II objects 




Class II 


277 


0.83 


43.41 


0.33 ±0.07 


17516 


SED II 


235 


1.43 


43.88 


0.12 ±0.07 


11305 



Bright23 

Bright23 23 1.13 44.77 0.84 ± 0.04 16281 
Note — (1) Samples of X-ray sources (see § 3 for details); 
(2) Number of sources; (3) Median redshift; (4) Median 2-10 
keV luminosity in logarithmic units of erg s~^; (5) Energy 
index of the stacked spectrum; (6) Number of net source 
counts. 




100 

Net counts 

Fig. 1. Distribution of the background-corrected counts of 
individual spectra in the rest-frame 2.1-10.1 keV band, ob- 
tained from XMM-Newton EPIC pn. 




Redshift 

Fig. 2. Distribution of redshifts for 1020 sources in the main 
sample (black). The Type I objects (blue) and Type II ob- 
jects (red) are also shown separately. The median redshifts 
for the total, Type I objects and Type II objects are 1.30, 
1.60, 0.97, respectively. 
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Fig. 3. Distribution of sources in the rest-frame 2-10 keV 
luminosity-redshift plane for a) Type I and b) Type II ob- 
jects. Sources in the main sample are indicated in black and 
the 23 sources in the bright sample are in red. Luminosity 
boundaries for investigating spectral variations are indi- 
cated by dashed lines (§ 5.1.1, § 5.2.1). Redshift bound- 
aries at z = 1 and z = 2 are also indicated for the Type II 
objects. 



Merloni et al 2010). Merloni et al (2010) use Mg ii while 
Trump et al (2009) give estimates based either on C iv, Mg 
II or Hz? where available. There are 37 objects for which 
both papers have measurements based on Mg ii. A com- 
parison between the two papers finds no systematic differ- 
ence in A/bh values, but there is a scatter of ~ 0.3 dex, 
which can be considered as the uncertainty in Mbh esti- 
mates. In the analysis of Trump et al (2009), A/bh based 
on Mg II is considered to be most robust, and in the event 
that a Mg ii based value is not available, estimates of A/bh 
based on other emission lines are adopted. Out of 366 Class 
I AGN in the main sample, 181 objects are found to have 
A/bh estimates. Given the availability of spectral lines, a 
large fraction of these objects (149) lie in the redshift range 
z =1-2.2. 

The accretion rate relative to the classical Eddington 
limit (against electron scattering), i.e. the Eddington ra- 
tio, A, is estimated from a ratio of the bolometric lu- 
minosity (Lboi) and the Eddington luminosity (£Edd — 
1.3 X 10^8 (A/bh /Mq) erg s'^). For almost all Class I AGN 
with determined A/bh (177 out of 181), their bolometric 
luminosities were computed by Lusso et al (2010) using the 
multi-wavelength SED. For these 177 objects, the range 
of log (A/bh/A/q) is between 7.18 and 9.70 with a me- 
dian value of 8.39 (Fig. 5a) while log (A) lies in the range 
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Fig. 4. Distribution of the rcst-framc 2-10 kcV luminosity 
of the sample sources: a) Type I, and b) Type II objects. 
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Fig. 5. Left: The distribution of the black hole masses es- 
timated for the spectroscopically confirmed Type I AGN 
(177 objects) taken from Trump et al (2009) and Merloni 
ct al (2010). The median value is 8.39. Right: The distri- 
bution of the Eddington fraction A for those Type I AGN. 
The median value is —1.0. 



— 1.9 to +0.4 with a median value of —1.0 (Fig. 5b). There 

is a weak inverse correlation between the accretion rate 
and black hole mass (see also Trump et al 2009; Netzer 
& Trakhtenbrot 2007), but with a large scatter, which is 
partly due to a selection bias in a flux-limited sample. This 
subset has a median redshift, z = 1.51, and median 2-10 
keV luminosity log (L2-10) = 44.26 crgs^^, and median 
bolometric luminosity log (^boi) = 45.52 ergs~^. 



b) 



c) 



d) 



6 6.5 7 
(keV) 



6 6.5 7 
(keV) 



6 6.5 7 
(keV) 



6 6.5 7 
(keV) 



Fig. 6. A narrow Gaussian was simulated at redshifts; a) 
z = 0, b) 0.5, c) 1 and d) 2, as observed with the EPIC pn, 
and then converted to a rest-frame spectrum after correct- 
ing for detector response and redshift for individual sources 
in the stacking sample. The rest-frame line energy is accu- 
rately recovered. The line profiles of higher redshift objects 
broaden because the energy resolution, (AE/E), of the de- 
tector becomes worse at lower energies. For this reason, 
while the integrated line flux is recovered, the line peak is 
lowered. 



3.3. Bright 23 

The 23 sources with more than 400 net counts each in the 
rest frame 2-10 keV band are mostly broad-line AGN (19 
Type I objects, 4 Type II objects) at ^; < 2. In Fig. 3, they 
are found in the upper envelope of the source distribution 
and represent the most luminous objects at a given redshift. 
The median properties are given in Table 1. 



4. Stacking method 

The rest-frame 2.1-10.1 keV data for individual sources are 
stacked. An original EPIC pn spectrum has 4096 channels 
with 5 eV channel width. The rest-frame spectral stack- 
ing is designed to have a 2.1-10.1 keV rest-frame energy 
range with 200 eV channel width. In an individual source 
with a redshift z, the EPIC pn data in the [2.1/(l+z)] keV 
to [10.1/(l-|-z)] keV range are selected and divided into 
40 energy intervals with [200/(l-|-z)] eV. These data, af- 
ter background subtraction, are corrected for the detector 
response curve and for Galactic absorption, A^h= 2.6 x 10^° 
cm~^(Kalbcrla et al 2006), before correcting the energy 
scale for redshift. 

The data stacking is a straight sum of individual 
sources. The reason why normalizing each spcxitrum by its 
flux was not made was the very large uncertainty in choos- 
ing the proper normalizing factor for faint sources, e.g., 
those with ~ 20 counts, which represent the majority of 
the sample. Thus the stacked spectra presented below nat- 
urally give more weight to brighter sources. However, our 
primary objective is to detect faint signals of Fe K emission 
in distant active galaxies that cannot reliably be detected 
by an ordinary observation of a single source, rather than 
to obtain a "mean" spectrum. 

As the energy resolution, AE/E, of the EPIC pn de- 
grades towards lower energies (approximately oc E~°-^), 
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a redshift-corrcictcid line profile of a higher redshift object 
would be broader than that of a lower redshift object. This 
effect is demonstrated in Fig. 6. Narrow Gaussians (with 
dispersion a = 1 cV) were simulated using an identical 
normalization, but at different redshifts, z =0, 0.5, 1, and 
2, and then the same correction procedure described above 
was applied to plot them in rest-frame energy. 

The energy resolution of the EPIC pn is « 0.15 keV 
(FWHM) at 6.4 kcV. However, given the design resolution 
of 0.2 keV, a simulated narrow line profile for z=0 - (a) 
has FWHM = 0.21 keV corresponding to a Gaussian dis- 
persion, a = 0.09 kcV. As redshift increases, the line profile 
broadens to a — 0.15 kcV at z = 2. While the Fe K line 
intensity distribution as a function of redshift is, of course, 
not known, if we assume that the median value of redshift 
is representative, the artificial line broadening due to in- 
strumental effects would be w 0.094 kcV for the stacked 
spectrum of Type I AGN, and a ~ 0.076 keV for Type 
II AGN. Thus, our integrated spectra are insensitive to a 
line width with ct < 0.1 keV in the Fe K band. The Fe 
lines discussed below are unresolved in most cases, other- 
wise the given resolution has been corrected for artificial 
broadening, assuming a typical redshift. 



5. Results 

For convenience of visual inspection, stacked spectra ana- 
lyzed in this paper arc presented in units of flux density 
{Fe), and the spectral slope refers to an energy index, a^, 
defined as Fe oc E~"'^. The photon index, F, commonly 
used in X-ray astronomy, is related to the energy index as 
F = 1 + QfE. Fig. 7 shows integrated (stacked) spectra for 4 
subsets of objects: (1) the main sample, (2) the "Bright 23", 
(3) Type I objects, and (4) Type II objects. 

The spectrum of the main sample has a spectral slope, 
aE = 0.60 ± 0.03 (hereafter errors quoted to spectral pa- 
rameters are of Ict) . This continuum represents a sum of all 
the sources in the rest-frame 2-10 keV band, and is steeper 
than the value of a-E ~ 0.4 for the X-ray background in 
the observed 2-10 keV band. A likely explanation for this 
difference is that bright, less obscured sources dominate 
the total integrated spectrum, given the sensitivity limit of 
the XMM-COSMOS. A significant Fe K feature is found at 
^ 6.5 kcV. Fitting the Fe K excess with a single Gaussian 
gives a line centroid of 6.55 ± 0.16 keV with a = 0.3lg'^ 
keV and an equivalent width, EW = 0.18±0.07 keV. When 
modelled with two Gaussians, which is more realistic, the 
line centroids are 6.36^° Qg keV and 6.86lg;j2 ^eV, with a 
flux ratio of 1:0.8. The former is unresolved and can be iden- 
tified with cold Fe K. The latter is marginally resolved with 
a = 0.15^0 A keV, and together with the centroid energy, 
this suggests a blend of Fe XXV at 6.70 keV and Fe xxvi at 
6.97 keV. No clear sign of a broad red wing is found. 

The integrated spectrmn of the "Bright23" sample 
shows a power-law spectrum with Oe = 0.84 ± 0.04. No 
signiflcant Fe K emission is detected. If anything, a broad 
emission feature with EW ~ 0.08 keV might be present, 
but its excess above the continuum is barely detected at 
the 1(7 level. The 2a upper limit for a narrow line at 6.4 
keV corresponds to EW = 0.09 keV. 

The spectra of Type I and Type II objects are discussed 
further in the following subsections. The spectrum of Type 
II objects is significantly harder {as — 0.2), as expected 




2 5 10 

Rest energy (keV) 



2 5 10 

Rest energy (keV) 



Fig. 7. Stacked spectra of the main sample, "Briglit23" 
sample. Type I objects and Type II objects. The spectra 
are presented in flux density units of 10"^^ [erg cm~^ s"""^ 
keV~^]. The rest energies of 6.4 keV and 6.9 keV are indi- 
cated by dotted lines. 



for various absorbed continua that are superposed, rather 
than that of Type I objects (aE — 0.8). 

To assess the effect errors in the photometric redshifts, 
the stacked spectrum for SED II sources is inspected, since 
these objects are expected to emit a narrow Fe K line. 
According to Salvato et al (2009), the accuracy of the pho- 
tometric redshifts is estimated to be ~ 1.5 %, which can 
be translated to tr ~ 0.1 keV for a narrow emission-line at 
6.4 keV. The stacked spectrum for SED II sources shows 
an unresolved (cr = O-ltH keV) line at 6.38 ± 0.13 keV 
with EW = 0.10 ± 0.07 keV. Relative to the spectroscopic 
sample of Type II objects (Class II), the EW of the 6.4 
keV emission is ^ 60%. Given the loose constraint on the 
line width, large errors in photometric redshifts cannot be 
ruled out as the cause of the reduction in the EW of Fe K. 
However, the best-fit line width agrees with the expected 
line broadening from the error reported in Salvato et al 
(2009), and the redshift/luminosity dependence (see § 5.2) 
may be a more likely cause, as the SED II sample includes 
more luminous, higher redshift objects than the Class II 
sample (see Table 1). 

5.1. Type I AGN 

Type I objects consist of 366 sources with spectroscopic 
redshifts (Class I) and 142 sources with photometric red- 
shifts (SED I). The sources in the photometric sample are, 
on average, fainter than those in the spectroscopic sample, 
and their contribution to the stacked spectrum is 20%. 

The 2-10 keV contimmm has a slope, aE = 0.80 ± 0.05, 
similar to the average value observed in nearby Seyfert 
1 galaxies and QSOs observed with XMM-Newton, e.g., 
the CAIXA sample (Bianchi et al 2009). The spectrum 
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of the photometric sample shows a harder continuum, 
a-E = 0.59 ± 0.07, than that of the spectroscopic sample, 
aE = 0.80 ±0.03. This suggests that the SED I sample con- 
tains some absorbed sources that might be mis-classificd as 
Type I AGN. The contamination rate was estimated to be 
12 % in Lusso et al (2010) for the objects for which both 
spectroscopic and photometric redshifts are available. 

There is excess emission at ~ 6.6 keV (Fig. 7), which 
appears to consist of two unresolved emission-line compo- 
nents at 6.4 ± 0.2 keV and 6.8 ± 0.2 keV rather than a 
single broad line. The line flux ratio of these two lines is 
approximately 2:3. The EW of the lines are 0.04 ± 0.02 
keV and 0.05 ± 0.02 keV, respectively. With a median 2-10 
keV luminosity of log L2-10 — 44.2 erg s~^, the EW of the 
6.4 keV line agrees with that expected from the Iwasawa- 
Taniguchi (IT) relation (Iwasawa & Taniguchi 1993) found 
for the CAIXA sample (Bianchi ct al 2007). The 6.8 kcV 
component is likely a blend of Fe XXV (6.70 kcV) and Fc 
XXVI (6.97 keV). These high-ionization lines are often ob- 
served in nearby AGN (e.g., Bianchi et al 2009) but they 
are usually weaker than the 6.4 keV line. This is reversed 
in our stacked spectrum of more powerful AGN at higher 
redshift.The nature of the enhanced high-ionization line is 
examined in § 5.1.2. 

No clear sign of redward broadening in the Fe K line 
profile, as observed in nearby Seyfert galaxies (e.g., Tanaka 
et al 1995), has been found either in the stacked data of 
the Type I objects, the Class I objects, or the SED I ob- 
jects. However, an upper limit for a relativistically broad- 
ened Fe K line with a shape typically found in nearby AGN 
(e.g., de la Calle Perez et al 2010), is ~ 0.2 keV. 

With regard to Fe K emission in Type I objects, we 
would like to highlight the following two points, which 
warrant further discussion since they may illustrate pos- 
sible differences in spectral characteristics between the dis- 
tant population of active galaxies, as sampled by XMM- 
COSMOS, and nearby active galaxies: 1) No strong line- 
broadening, such as might result from relativistic effects is 
found; and 2) The high-ionization lines Fe XXV and Fe xxvi 
increase in strength. 

In order to search for a population that emits strong 
broad iron lines, or the high-ionization lines, we explore the 
Type I dataset by dividing it into three subsets according 
to the 2-10 keV luminosity, redshift, and black hole mass 
or Eddington ratio. Here we report our findings in terms of 
the above two points. 

5.1.1. Broad Fe K line 

To avoid a possible artifact due to errors in redshift esti- 
mate, which could result in a false detection of line broad- 
ening, we discarded objects with photometric redshifts and 
used only the spectroscopic sample (Class I) for the further 
investigation of broad emission. Empirically, it has been 
found that at least ~150,000 counts in the 2-10 keV band 
are required for a significant detection of a broad Fe K line 
with EW ~ 0.1-0.15 keV, typical of that found in nearby 
Seyfert galaxies, with the EPIC camera (dc la Calle Perez et 
al 2010). The integrated counts of the entire sample of Class 
I objects is ~40,000, well below the critical counts (the 
background is also much higher), and only strong broad 
emission, with EW >0.2-0.3 keV, can be detected. Bearing 
this limitation in mind, nevertheless, it is found that there 
may be a favourable range in X-ray luminosity and/or in 
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Fig. 8. Stacked spectra of Class I objects (BLAGN) in the 
luminosity ranges above and below log £2-10 = 44.43 (erg 
s~^). For the spectrum of the high luminosity objects, the 
best-fit power-law continuum is indicated by the solid-line 
histogram. Weak Fe K emission peaking at 6.86 keV is seen. 
On the other hand, the continuum model indicated with the 
solid-line histogram for the spectrum of low luminosity ob- 
jects is that expected for a power-law with reflection (see 
text for detail). A moderately broadened Fe K feature cen- 
tred at 6.5 keV, with a possible redward extension down 
to 4.5 kcV, can be seen. The best-fitting diskline model for 
this Fe K feature is indicated by the dotted-line histogram. 

Eddington ratio for finding broad Fe K emission, the details 
of which are described below. 

When the Class I objects are divided into two X-ray 
luminosity ranges, above and below 3 x 10^** erg s~^, a 
possible signal of a broadened Fe K emission line is found 
in the lower X-ray luminosity subset, albeit the significance 
of the broad component is barely at 2(t level. This subset 
consists of 247 objects with median redshift z — 1.45 and 
median 2-10 keV luminosity, log (^2-10) — 44.00 erg . 

The stacked spectrum is shown in Fig. 8. The continuum 
slope is Q!E = 0.88 ± 0.05. Fitting a single Gaussian to the 
Fe K feature gives a line centroid of 6.55 ± 0.26 kcV with 
a = 0.37to.i6 keV. This simple fit leaves a weak but positive 
residual down to 4.5 keV, suggesting the presence of a red 
wing. For comparison, the spectrum of objects in the higher 
luminosity range is also shown. The spectrum of the high 
luminosity subset agrees well with a simple power-law of 
cce = 0.80 ± 0.05 apart from narrow emission-line features 
at 6.86 ± 0.15 keV and possibly at 6.38 ± 0.21. There is 
no hint of excess emission below 6 keV. This subset carries 
a factor of 1.4 more flux than the lower luminosity subset 
with the broad feature, which would explain why no broad 
emission was seen in the total spectrum. 

While the line broadening is marginally significant, pres- 
ence of a red wing below 6 keV is highly uncertain and 
introducing any further complexity in a spectral fit will be 
subject to overmodelling. However, in the interest of a com- 
parison with the nearby AGN, we present a modelling of the 
broad feature, assuming that it is an emission line arising 
from a relativistic disk. A reflection continuum could mimic 
a weak red-wing and its effect appears at maximum without 
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rclativistic smearing, as modelled with pexrav (Magdziarz 
& Zdziarski 1995). The broad feature is found to be robust 
aginst this. In the continuum model, the cold reflection is 
included with a reflection fraction of unity and inclination 
of the reflecting surface to be 35° (see below) . The line fea- 
ture is modelled by the relativistic line model diskline 
(Fabian et al 1989). Given the quality of our data and the 
degenerate nature of some model parameters, typical val- 
ues for the disk-line model were chosen and only the disk 
inlination and the line normalization are fitted. The disk 
radii were assumed to be between 10-20 Vg (see Nandra et 
al 2007 for the typical radius) with an emissivity index of 
—2, and the line was assumed to be emitted at 6.4 keV. 
Combined with the continuum model including reflection, 
the best-fit inclination is found to be 35° ± 7° . The line flux 
corresponds to EW = 0.28l^:i2 keV. This model repro- 
duces the moderately broad 6.55 keV peak, picked up by 
the single Gaussian fit and the possible redward extension 
down to 4.5 keV. No rapid spin of the black holes is inferred 
from the redwards extension of the line in the present data, 
but much higher quality of data in the red-wing band are 
needed to impose a meaningfiil limit on the spin. 

The above EW is likely an upper limit for any broad 
component of the line, since the emission peak at 6.55 keV 
is at least partly due to a blend of the cold Fe line at 6.4 keV 
and high-ionization lines at 6.7-7 keV. It has been known 
that a narrow Fe K line at 6.4 keV, originating from cold 
matter at large distances, is ubiquitously found in AGN. 
Using the median X-ray luminosity of the lower luminosity 
subset, combined with the known IT relation (Bianchi et al 
2007), the EW of the 6.4 keV line is expected to be ~ 0.05 
keV. Furthermore, a contribution by the high-ionization 
lines with EW ~ 0.05 keV is also expected at energies 
corresponding to the blue wing. Thus, if present, the true 
contribution of the broad component would be EW < 0.2 
keV. Note also that the subtraction of the non-relativistic 
components affects the diskline parameter, in particular the 
inclination, which would be lowered. 

A similar broad emission feature can also be found when 
an intermediate range of the Eddington ratio, e.g. A ^ 0.1, 
is selected (Fig. 9). The details of the source selection ac- 
cording to A is presented in Table 2. Although the signal 
to noise ratio is lower, given that the number of objects in 
the stack is only 97, besides a narrow line at 6.50 ± 0.02 
keV, a red wing can be seen down to 5 keV. The level of 
detection of this broad feature is also ~ 2a. The typical X- 
ray luminosity (log ^2-10 — 44.17 ergs~^) of this subset is 
within the low X-ray luminosity range (log 1/2- 10 < 44.43 
erg s^^) where the broad emission seems to be present. This 
makes it difficult to determine which parameter drives the 
production of possible broad Fe K emission. 

In summary, a broad red wing of Fe K emission in Type 

I AGN may be present, when the X-ray luminosity range 
is restricted to below ~ 3 x lO''^ ergs^^, or when sources 
with an intermediate range of Eddington ratio (—1.3 < 
log (A) < —0.6) are selected. The detection of the broad 
emission is, however of low significance, and the EW of the 
broad feature would be < 0.2 keV, in agreement with the 
previous studies by Corral et al (2008) and Chaudhary et 
al (2011). Note that our objects lie, on average, at higher 
redshifts than those in the previous work, and the sample 
of Chaudhary et al (2011) contains both Type I and Type 

II AGN, which makes a direct comparison difficult but an 
overall agreement on broad emission is found. As argued 
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Fig. 9. Possible broadened Fe line profiles found in spectra 
of a) lower luminosity range in Class I objects and of b) the 
intermediate range of Eddington ratio (Table 2 and Fig. 
10b), presented in the form of ratio versus the power-law 
contimmm. The data are rebinned to 600 eV intervals for 
display purposes. 



above, our data are expected to be only sensitive to strong 
broad-lines with EW > 0.2-0.3 keV, and the obtained limit 
is as expected. 

5.1.2. High-ionization Fe emission 

An investigation of the Fe K line properties using the Class 
I objects suggests that the Eddington ratio, A, may play a 

role for producing the high-ionization Fe K, which is seen 
strongly in the spectrum of Type I objects (§ 5.1 and Fig. 
7). 

In the following discussion we divide the 177 objects 
with A estimates (§3.1.1, Fig. 5b), into three subsets accord- 
ing to A: 1) Al (log a < -1.3), 2) Am (-1.3 < logA < -0.6), 
and 3) Ah (log A > —0.6). The properties of the three sub- 
sets are given in Table 2. Fig. 10 shows the stacked spectral 
data from the three A intervals. The stacked spectrum of the 
objects in the high A range (Ah) shows an iron line at 6.9 
keV. Prior inspection of individual spectra of the brightest 
sources in this subset shows that a single source, XID 5230 
might dominate the line feature. The EPIC pn spectrum of 
this source, as observed, shows a clear line feature at an ob- 
served energy of 3.00±0.03 keV, which corresponds to 6.95 
keV in the rest-frame of this galaxy {z = 1.32). Therefore 
this source was excluded from the stack (i.e. the spectrum 
in Fig. 10 was constructed excluding XID 5230). However, 
despite the exclusion of XID 5230, the stacked spectrum of 
the remainder of the sources still shows a line feature at 
6.89^^:^^ keV with EW 0.19 ± 0.10 keV. The fine flux is 
comparable (~75%) to that solely from XID 5230. A cold 
Fe K line at 6.4 keV is not detected, with a 2<t upper limit 
of EW = 0.07 kcV. 

Besides the high accretion rate, the objects with Ah, 
exhibiting the high-ionization Fe K line, are characterized 
by greater X-ray output (log L2-10 = 44.52 erg s"-*^) from 
relatively low-mass black holes [log {Mqu/Mq) = 8.1] re- 
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Table 2. Properties of Class I objects in three ranges of accretion rate, A. 



Subset 


Range 


N 


Mbh 


A 


z 


^2-10 


F 


F/N 


Oe 


Fe K 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Al 


< -1.3 


40 


8.51 


-1.48 


1.20 


44.18 


4.2 


1.1 


0.87 ±0.10 




Am 


-1.3 - -0.6 


97 


8.39 


-1.00 


1.51 


44.17 


9.3 


0.96 


0.86 ± 0.07 


Broad? 


Ah 


> -0.6 


39 


8.10 


-0.31 


1.82 


44.52 


4.5 


1.2 


0.77 ± 0.07" 


6.9 keV 



Note — 177 objects with lambda estimates are divided into three subsets according to their value of A. One source XID 5230 in 

the Ah subset is excluded (see text). 

(1) Subset according to the classical Eddington ratio, A. 

(2) A range in logarithmic unit. 

(3) Number of sources. 

(4) Median black hole mass in logarithmic solar mass unit. 

(5) Median Eddington ratio A in logarithmic unit. 

(6) Median redshift of the sources. 

(7) Median rest-frame 2-10 keV luminosity of the sources in logarithmic units of erg s~^. 

(8) Integrated rest-frame 2-10 keV flux of the stacked spectrum in logarithmic units of 10~^^ erg cm~^ s"'^. 

(9) Arithmetic mean of the 2-10 keV flux per source in logarithmic units of 10^^'' erg cm^^ s^^. 

(10) Spectral slope of the stacked spectrum in the rest-frame 2-10 keV band, a: When absorption is included, a steeper value 
QE = 0.99 ± 0.19 is obtained with A'h = (0.85 ± 0.64) x 10^^ cm"^. 

(11) Comment on the Fe K feature. 
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Fig. 10. Stacked spectra of BLAGN (Class 1 objects) in three A(= iboi/^Edd) ranges: a) log A < —1.3; b) —1.3 < log A < 
—0.6 and c) log A > —0.6 (see Table 2 for details). The rest energies of 6.4 keV and 6.9 keV are indicated by dotted lines. 
The stacked spectrum of objects with high accretion rates shows a strong Fe K line at 6.9 keV. 



siding in galaxies at slightly higher redshift (z = 1.82), rel- 
ative to those with lower values of A. These characteristics 
are also shared by the excluded source XID 5230, which 
has Mbh = 8.21, log A = -0.04, log L2-10 = 45.68 and 
z = 1.32. Note that 15 sources out of 39 in the Ah subset 
have Mbh < IO^Mq. This implies that, when a stacked 
spectrum is made for objects with log (Mbh/M©) < 8, the 
spectrum also has a 6.9 keV line. 

While no significant line is detected, the presence of 
the 6.9 keV line cannot be ruled out for the noisy Al 
stacked spectrum. The EW of the 2fT upper limit is 0.35 
keV. Likewise, any broad emission would have EW < 0.6 
keV. A tight upper limit for a 6.9 keV fine with EW < 0.08 
keV, is obtained for the Am spectrum. 

There is a weak trend corresponding to an anti- 
correlation between Mbh and A, as noted in § 5.1.1. To 
verify that A, rather than Mbh, indeed drives the spec- 
tral variations seen above, we repeated the same exercise 
in a common range of Mbh- Limiting the Mbh range to 
1O^-°M0 - 10^-^ Mq, i.e., ±0.4 dex around the median value 
(see Fig. 5), the numbers of objects contained in the three 
A ranges change to 31 (Al), 77 (Am), and 21 (Ah, XID 5230 
excluded). The median log (A'/bh) is now found to be 8.4 
for all the three A intervals. Spectral properties obtained 



from these data agree well with those without the Mbh 
restriction, except for a slight steepening of the continuum 
slope, Q!E = 0.95±0.09, and an increased EW = 0.27±0.14 
keV for the Ah subset. This verifies that A is indeed a likely 
driver of the variation in the Fe line characteristics found 
above. 



5.1.3. Spectral slope 

Previous studies have shown that the X-ray spectral slope 
steepens with increasing Eddington ratio, A (e.g., Shemmer 
et al 2006; Porquet et al 2004). This trend is not clearly 
seen in our Class I objects (Table 2). The spectrum for 
the Ah objects might have a curvature at lower energies 
which could be attributed to extra absorption. If absorp- 
tion is included, the continuum slope would steepen to 
aE = 0.99 ±0.19, with A^h = (0.9 ±0.6) x 10^2 cm'^. This 
slope is similar to that (~ 0.95) obtained when the Mbh 
restriction is applied (without extra absorption). Even tak- 
ing this steepened slope, a change in slope AaE ~ 0.1 over 
the 1.2 dex interval in A between the medians of Al and 
Ah appears smaller than that previously reported, e.g., in 
Shemmer et al (2006). Conversely, if an A relation is 
true, then the apparently similar slopes between the three 
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Fig. 12. Plot of EW of the Fe K line at 6.4 keV (open 
circles) and 6.9 keV (filled squares) measured from a) the 
three 2-10 keV luminosity bins or b) the three redshift bins 
(see Table 3). The EW is plotted against the median value 
of luminosity/z in each bin. Note that the high-z point for 
the 6.9 keV line is a sum of Fe XXV (6.70 keV) and Fe xxvi 
(6.97 keV). 

subsets would suggest increasing absorption with increas- 
ing A, which would cancel out any such steepening of the 
slopaM Such absorption may be due to a partially ionized 
(warm) medium so that the visibility of the optical BLRs 
are not affected. 

5.2. Type II objects 

The stacked spectrum of Type II objects shows a hard con- 
tinuum (a = 0.22 ± 0.05) with a strong Fe K feature (Fig. 
7). The Fe K feature is resolved into two components at 
6.42j:o;o7 keV and at 6.9ltom keV, with a line intensity 
ratio of 1.0 : 0.4. The two lines are both unresolved individ- 
ually and have EW = 0.15±0.05 keV and EW = 0.08±0.06 
keV, respectively. We refer to Mainieri et al (2011) who 
present the X-ray results and the host galaxy properties 
for a subset of this sample, i.e. Type II QSOs with high 
luminosity and significant X-ray absorption. 

The spectral properties are also investigated by divid- 
ing the Type II dataset into three ranges of the 2-10 keV 
luminosity and three ranges of redshift (Fig. 11); their ba- 
sic properties are listed in Table 3. The continuum slope 
remains nearly constant, a ~ 0.25 ± 0.07, except for the 
spectrum of the ZH {z > 2) sample where a = 0.14 ± 0.11. 
Interesting variations in Fe K emission were also found, as 
reported below (see Fig. 12). 

5.2.1. Luminosity dependence 

The three luminosity intervals (LL, LM and LH) are shown 
in Table 3 (see also Fig. 3b) and the stacked spectra of 
sources in those ranges are shown in Fig. 11a, b,c, respec- 
tively. The Fe K feature weakens as the 2-10 keV luminosity 
increases (Fig. 12a). The cold line at 6.4 keV (open circles 

^ A spectral curvature due to absorption appears in a stacked 
spectrum only when the absorbing columns are similar between 
objects, otherwise they effectively make a power-law slope flat- 
ter. 



in Fig. 12a) is seen in all three luminosity ranges but its EW 
decreases with increasing luminosity, similar to what is seen 
in Type I AGN. The high-ionization line at 6.9 keV shows 
a similar decreasing trend (filled squares in Fig. 12a), with 
a line flux ratio of ~ 0.7 compared to the cold line. This 
trend can be compared with the earlier result by Krumpe 
et al (2008). 

5.2.2. Redshift dependence 

Stacked spectra for sources in the three redshift intervals, 
ZL (z < 1), ZM (1 < z < 2) and ZH (z > 2) (see Table 3 
and Fig. 3b) are shown in Fig. lld,e,f, respectively. Only 
the spectrum for objects with z < 1 shows the 6.4 keV 
line as well as the high-ionization line at 6.85 keV. In the 
intermediate interval (1 < z < 2), no significant line is 
detected, although a weak excess is seen around 6.5 keV 
(EW < 0.15 keV). However, at z > 2, a significant Fe K 
feature reappears. There are sufficient statistics to allow 
this feature to be resolved into two components: Fe xxv 
(6.70 keV) and Fe xxvi (6.97 keV) with EW = 0.21 ± 0.1 
keV and EW = 0.13 ± 0.07 keV, respectively. No 6.4 keV 
line is seen (the 2a upper limit for the EW is 0.1 keV). 

One interpretation for the absence of a line feature in 
the z = 1-2 stacked spectrum is that the iron atoms are 
predominantly in intermediate ionization states where the 
fluorescence yield is low (due to resonant trapping, e.g., 
Krolik & Kallman 1987), before moving up to higher ion- 
ization stages, as seen at higher z. Alternatively, since a 
large number of objects in this redshift range have only 
photometric redshifts (127 out of 185) any emission feature 
might be blurred. However, the z > 2 subset (ZH) also has 
a high fraction of objects with photometric redshifts (47 
out of 60), but the Fe line is clearly detected. 

As a remark related to § 5.2.1, although z and L2-10 
are, in principle, correlated, the zones occupied by the high 
-^2-10 (LH in Table 3) and the high z subsets (ZH) in the 
z-ix plane substantially differ from each other (see Fig. 3b 
and Table 3); their median luminosities are similar but the 
median redshifts are 1.76 (LH) and 2.52 (ZH), thus ZH is 
practically the high-z subset of LH. We discuss the origin 
of the clear appearance of the high-ionization Fe line when 
the redshift range is moved to z > 2 (§6.2 and §6.3). 

6. Discussion 

6.1. Broad Fe emission 

Broad Fe K emission has been found in spectra of nearby 
unobscured AGN. Given that the primary X-ray production 
likely occurs near the innermost radii of the accretion disk, 
and illumination of the disk causes the line emission, rela- 
tivistic broadening is considered to be a likely explanation 
(e.g., Fabian et al 2000), while a complex absorption model 
has been proposed to explain the redward asymmetry (e.g., 
Miller, Turner & Reeves 2009). Whatever the origin of this 
broad feature, previous investigations of Fe line profiles of 
nearby, bright Seyfert galaxies found that the fraction of 
objects exhibiting such broadening is /b ~ 40 %, for which 
the mean EW of the broad lines is EW ~ 0.1 — 0.15 kcV 
(Nandra et al 2007; de la Calle Perez et al 2010). This im- 
plies that the mean spectrum of these galaxies would have 
a broad line with EW = 0.04 - 0.06 keV. Our adopted 
stacking method means that a direct comparison cannot 



10 



K. Iwasawa et al.: Fe K emission of active galaxies in COSMOS 




CM 



d 



b) LM 






1 : 

Iff 



U3 

d 



10 



e) ZM 


.1 1 1 1 




■, , , 1 



in 



CM 





.1 1 1 1 

: 


>++^%tt,tt|t||j 






;, , , r 



10 




2 5 10 

Rest energy (keV) 
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Table 3. 
jects. 



Luminosity and redshift intervals of Type II ob- 



Subset 


Range 


L 


z 


N 


F 


FIN 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




log L2-10 












LL 


41.00-43.35 


43.05 


0.61 


162 


0.40 


0.25 


LM 


43.35-43.95 


43.68 


1.03 


211 


0.76 


0.36 


LH 


43.95-45.50 


44.22 


1.76 


139 


1.30 


0.94 


ZL 


z 

0-1 


43.31 


0.71 


267 


1.14 


0.43 


ZM 


1-2 


43.84 


1.40 


185 


0.91 


0.49 


ZH 


2-5 


44.30 


2.52 


60 


0.43 


0.71 



(1) 

(2) 
(3) 
(4) 
(5) 
(6) 
ergs ^ cm 

(7): Mean rest-frame 2-10 keV flux per source in unit of 



Luminosity or redshift subset. 

Range of variable, either 2-10 keV luminosity or redshift. 
Median 2-10 keV luminosity. 
Median redshift. 

Number of sources contained in the subset. 
Integrated rest-frame 2-10 keV flux in unit of lO"^'^ 



10" 



ergs 



be made, but this value is far below our detection limit. 
Therefore, no detection of broad emission is perhaps not a 
surprise. 

On the other hand, when the low luminosity range or 
intermediate Eddington ratio are selected, possible broad 
line signals are seen (sect. 5.1.1). If the broad feature with 
EW ^ 0.2 kcV is real, strong, broad Fe emission could be 
ubiquitous (/i-, ^ 1) among Type I AGN with ^2-10 ~ 10*^ 
ergs^^ and A ^ 0.1 at z ^ 1.5. This has to be tested 
by observations of individual galaxies that meet those se- 



lection characteristics, and may give impetus for a future 
large X-ray telescope mission. It is also in line with the 
theoretical prediction made by Ballantyne (2010), although 
the favourable luminosity range is lower. In nearby AGN, 
objects with high A like MCG-6-30-15 (e.g., Tanaka et al 
1995) seem to show broad Fe K emission more often than 
normal Seyfert 1 galaxies, which is contrary to our finding. 
However, in a nearby AGN sample of FERO (de la Calle 
Perez et al 2010), no A dependence in the detectabihty of 
broad emission is reported. 

6.2. High-ionization Fe K in distant AGN 

High-ionization Fe K emission from Fe XXV -|- Fe XXVI ap- 
pears to increase in importance compared to the cold line 
found at 6.4 keV in the spectra of the XMM-C0SM:0S 
AGN, compared to well-studied, nearby active galaxies. 
The ionized gas emitting these lines is in such a high ion- 
ization state that other lighter elements are expected to be 
almost fully stripped. These high-ionization lines are found 
to be the major component of the Fe K complex in the spec- 
tra of Type I AGN with high Eddington ratio (A > 0.25, 
§ 5.1.2), and Type II AGN at high redshift (z > 2, § 5.2.2). 
Here we discuss the origin of the high-ionization Fe K emis- 
sion in the two populations as well as a possible connection 
between the two. 



6.2.1. BLAGN with high Eddington ratio 

The investigation of BLAGN in § 5.1.2 suggests that there 
may be a causal link between the high-ionization Fe K and 
A. It has been reported that PG quasars with a high A at low 
redshift, e.g., z < 0.4, the same range as the CAIXA sam- 
ple, show high-ionization Fe K (Nandra et al 1996; Porquet 
et al 2004; Inoue, Terashima & Ho 2007) as well as some 
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narrow-line Seyfert 1 galaxies (Comastri et al 1998) that 
are also considered to have high A. 

With a high accretion rate, strong radiation heats the 
surface of the accretion disk and probably causes the inner 
part of the disk to inflate. This region is then expected to be 
highly ionized. The observed 6.9 keV line is narrow however, 
which excludes this relativistic region as the origin. Even in 
the absence of relativistic broadening, Fe emission produced 
in a highly ionized disk is broadened by Compton scatter- 
ings within the hot layer of the disk (e.g., Ross & Fabian 
2005). Further broadening by relativistic effects might make 
such line emission difficult to detect. 

Instead, the tenuous disk atmosphere which is likely 
outflowing due to X-ray heating (e.g.. Shields et al 1986) 
or torus winds which evaporate off the illuminated inner 
surface of the obscuring torus (Krolik & Begelman 1986; 
Balsara & Krolik 1993), could be a likely source of the high- 
ionization line. Such line emission is well predicted in both 
cases (Shields et al 1986; Krolik & Kallman 1987). The lack 
of a 6.4 keV line suggests that the covering factor of thick, 
cold clouds in these objects is small, although the 2(t upper 
limit is still consistent with that expected from the EW-Lx 
relation. 

6.2.2. High-z Type II objects: the local analogue and 
Eddington ratio 

The Type II objects at z > 2 (ZH) show an Fe K complex 
dominated by the high-ionization lines while a cold line at 
6.4 keV is undetected. We note that the hard continuum 
spectrum (a ~ 0.1), which is comparable or possibly even 
harder than those for the spectra in the lower luminos- 
ity/redshift ranges, indicates that a significant fraction of 
these ZH sources are strongly absorbed, yet the observed X- 
ray luminosity is very large (the median log ^2-10 — 44.30 
ergs""'^; see Table 3). This X-ray luminosity is only 2/3 
of the median value for the other high-ionization Fe line 
emitters among the Class I objects with high A. Combined 
with the EW, these two subsets are found to have super- 
ficially comparable line luminosities for the high-ionization 
Fe. This may indicate a common origin for the line emis- 
sion, and thus the high-z Type II objects contain black holes 
accreting at a high rate (or a high Eddington ratio. A). 

The quasar- like X-ray luminosity (> lO'*^ ergs~^), the 
hard continuum and the high-ionization Fe K of the high-z 
Type II objects are rarely seen in active galaxies in the local 
Universe. One of the few exceptions that have been known 
to have all of these characteristics is the highly obscured 
ULIRG, IRAS F00183-7111 {Li, ~ 7x at z = 0.33). 

It contains a heavily obscured AGN with Lx ~ 10**^ ergs~^, 
and shows a reflection-dominated X-ray spectrum with a 
strong Fe xxv line (Nandra & Iwasawa 2007). A remark- 
able characteristic of this ULIRG is the presence of fast out- 
flows (up to ~ 3000 km s^^) probed by optical [Oiii]A5007 
(Heckman, Armus & Miley 1990), and the mid-IR [Ne 11] 
lines (Spoon et al 2009). A small-scale radio jet imaged with 
VLBI supports the idea of an AGN driven outflow (Norris 
et al 2011). Although no black hole mass estimate (hence 
A) is available, a possible connection between the nuclear 
outflow and the production of high-ionization Fe K is sug- 
gested (e.g., Nandra & Iwasawa 2007). In local ULIRGs, 
molecular outflows are common and AGN radiation pres- 
sure seems to be a favoured driving mechanism (Feruglio et 
al 2010; Sturm et al 2011; but see Chung et al 2011). 



Lusso et al (2011) arc able to provide robust estimates of 
the bolometric luminosity estimates and stellar mass (Af^) 
for 267 Type II objects, based on their SEDs. Bearing in 
mind that M^, is not for a spheroid, but for an entire galaxy, 
we caution that a blind use of the scaling relation between 
M^, and Mbh is dangerous (resulting in an over-estimate 
of Mbh in most cases). Therefore, we attempted to make 
a crude estimate for A by using the relation obtained by 
Haring & Rix (2004). The median value for all 267 objects 
is A = 0.037. Unfortunately, only seven objects are found 
in the z > 2 range, but they all lie at the higher end of 
the distribution with a median value, A = 0.68 (practically, 
they are objects with higher Lx/Ai*). While the caveat 
above, the small number of the subsample, and the possible 
evolution of the scaling relation (Merloni et al 2010; Li, Ho 
& Wang 2011) make any results inconclusive, this indicates 
that at least some of these objects could indeed be obscured 
AGN accreting at a high Eddington ratio. 

6.3. HLIRGs at z>2 

Out of the 60 Type II objects in the ZH sample discussed 
above, mid-IR counterparts detected in the MIPS-24/xm S- 
COSMOS survey (Sanders et al 2007j^have been identifled 
for 27 objects (Brusa et al 2010). With a typical redshift 
oi z = 2.5 (Table 3), the detected emission corresponds to 
rest-frame emission at 7/zm with a median luminosity of 
log Li^rn = 45.20 (ergs~^). This emission probably comes 
from hot dust heated by an embedded AGN. When the 
mean SED of nearby ULIRGs derived from the GOALS 
sample (U et al 2011) is assumed, log Li, ~ 46.3 ergs~^, or 
^ 1O^'^L0, is obtained. 

These powerful mid-IR sources are found to be a dom- 
inant source of the high-ionization Fe K in the spectrum 
of the ZH sources (see Fig. llf). The stacked spectrum of 
these 27 MIPS-24/xm detected sources shows a Fe K com- 
plex with comparable intensity of Fe xxvi and Fe xxv (Fig. 
13). The Fe xxvi line accounts for virtually all the Fe xxvi 
flux in the stacked ZH spectrum. The spectrum of the other 
33 sources appears to only make a small contribution to Fe 
xxv (see Fig. 13). 

We note that several nearby (U)LIRGs, including Arp 
220, with no apparent AGN signatures, show strong Fe xxv 
(Iwasawa et al 2005, 2009), which can be accounted for by 
thermal emission originating from a starburst, as the X- 
ray to infrared luminosity ratio for those (U)LIRGs is log 
(Lx/iir) ^ —4.5 (Iwasawa et al 2009). In contrast, the 
ZH objects have log {Lx/Li,) ^ -2 and log Lx ^ 44.30 
(ergs~^). This X-ray loudness indicates that an AGN is a 
most likely origin for the X-ray luminosity and the high- 
ionization Fe K line in the stacked ZH spectrum. 

The number of known hyper-luminous infrared galaxies 
(HLIRGs: Ljr ^ 10^^ Lq) at 2: > 2 have increased dra- 
matically after the launch of the Spitzer Space Telescope. 
Many of the highly obscured HLIRGs at z > 2, selected at 
MIPS-24^m by Houck et al (2005), were found to have IR 
SEDs similar to that of F00183-7111 mentioned in § 6.2.2., 
which is characterized by little PAH emission and deep sil- 
icate absorption, suggesting a deeply buried AGN (Spoon 
et al 2007) . Four galaxies with their infrared SEDs that re- 
semble that of IRAS F00183-7111, were followed up by the 

^ The cataloK is publicly available thro ugh IRSA, 
http://irsa.ipac.caltech.edu/Missions/cosmos.html 



12 



K. Iwasawa et al.: Fe K emission of active galaxies in COSMOS 



m IT) 




o - 

_l I I I I I I I I I I I I I I I 

6 7 8 

Rest energy (keV) 

Fig. 13. The Fe K band stacked spectrum for 27 MIPS 
(24 /^ni)-detected objects (black solid squares) and 33 non- 
detected objects (red open circles) from the ZH subset 
(Table 3, Type II objects at z > 2). The stacked spec- 
trum of the MIPS detected objects shows much stronger 
line emission and these objects are a dominant Fe K source 
in the ZH stacked spectrum of Fig. llf. In particular, Fe 
XXVI appears to come only from the MlPS-detected objects. 



Chandra X-ray Observatory. Three of these sources were 
detected, with L2-10 ^ 10**^ ergs~^(Vignali et al, in prep). 
At z ~ 2, heavily obscured AGN are found to be abun- 
dant in BzK-selected galaxies with MIPS-24/im detections 
using an X-ray image stacking analysis (Daddi et al 2007; 
Alexander et al 2011), although they are a much less lumi- 
nous (1-2 orders of magnitude less) population. 

The MIPS-24/im detected Type II AGN at z > 2 are, 
on average, powerful infrared galaxies containing obscured 
AGN, as implied by the hard continuum of the stacked spec- 
trum (as — 0.1). If the connection between high- ionization 
Fe K and high A found in BLAGN (see § 6.2.1) is, also appli- 
cable here, it would mean that black holes in these objects 
are accreting at a substantial rate. The redshift interval 
where they lie coincides with the peak period of building 
massive galaxies and perhaps concurrent black hole growth 
(e.g. Chapman et al 2005; Hopkins et al 2006; Veilleux et 
al 2009). 

6.4. Compton thick AGN population at z ^ 0.7 

The Type II objects in the low luminosity (LL) and the low 
redshift (ZL) intervals have typical redshifts, z = 0.6 — 0.7, 
where the major contribution to the X-ray background is 
expected. Also, the Fe K features are found to be strongest 
in their spectra. A strong [EW ~ 1 keV), cold Fe line at 
6.4 keV is a characteristic signature of heavily obscured 
AGN. This large EW results from the suppression of the 
primary continuum due to strong absorption. Therefore, in 
our stacking method of straight integration, the Fe line flux 
comes from all sources while the continuum is dominated 
by less absorbed sources. The EW of Fe K in the LL and ZL 
stacked spectra is EW = 0.3 - 0.4 keV (see Fig. 12). This 
is slightly larger than the typical value {EW ~ 0.2 keV) 
of absorbed AGN with 7Vh= lO^^-lO^^ cm'^ (Fukazawa et 
al 2011), suggesting that there is some contribution from 



Compton thick AGN to the Fe K line. The LL and ZL 
sources have a typical flux level of 10^^"* ergs^^cm^^, 
about 2 orders of magnitude below that of bright Compton 
thick AGN detected with INTEGRAL IBIS or Swift BAT 
in the observed 2-10 keV band (e.g., Comastri et al 2010). 
Down to this flux level, no drastic change in the observed 
fraction of Compton thick AGN is expected and it would 
remain at 10-15 % (GiUi et al 2007; Treister, Urry & Virani 
2009 and references therein). The observed EW (Fe K) is 
consistent with this Compton thick AGN fraction. 
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